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1. Introduction 

There is observational evidence that the Universe is undergoing an accelerated 
expansion driven by a cosmological constant or some form of energy that violates the 
strong energy condition (p + 3p > 0, p and p being the energy density and pressure, 
respectively, of the fluid conforming the Universe) [TJ. The observational data seem 
to admit models that even assume that the Universe is dominated by a kind of fluid 
that violates the dominant energy condition (p + p > 0). This fluid is known in the 
literature as the phantom energy [2]. 

Many types of models have been proposed in recent years to explain the 
observed accelerated expansion: the cosmological constant, which presents a fine- 
tuning problem [3]; quintessence models [4] and quintessential inflation fields [5]; 
scalar-field models [6]; chameleon fields [7]; K-essence models [8 a ; modified gravity 
models [9] ; the feedback of non-linearities into the evolution equations [10] ; Chaplygin 
gases [IT] : tachyons [12]; phantom dark energy and ghost condensates [2]; de-Sitter 
vacua with the flux compactification in string theory |13] ; cyclic universe [14j . A more 
complete list and discussion can be found in Ref. [15] . Although all the models predict 
an accelerated expansion that is in good agreement with the observational data, few 
of them propose a microscopic understanding of dark energy. 

Another stage of accelerated expansion has been proposed over the years in 
an attempt to explain the horizon, flatness and unwanted relic problems: inflation 
[TBI [T71 ITS] . Most inflationary models assume that the expansion of the universe 
is dominated by a scalar field with a slow-roll potential, so that the scale factor 
evolves nearly exponentially with time. During the inflationary stage of the expansion 
the scale factor increases several orders of magnitude and the relic radiation energy 
density decreases dramatically. To allow nucleosynthesis to happen in the inflationary 
scenario, a reheating process should be introduced: the inflaton field decays into 
conventional fields by means of the potential. After reheating, the universe is 
dominated by radiation. 

Given the lack of information about the nature of the dark energy, 
thermodynamics and holography may be helpful. As examples of this approach, 
Pollock and Singh in Ref. [TjJ[ study the thermodynamics of the de Sitter spacetime 
and quasi-de Sitter spacetime. One of the quasi-de Sitter models studied presents 
phantom-like behavior and the entropy associated to it is negative defined. In Ref. 
[20] , an extension of these results is found for a more general phantom- type spacetime 
leading to negative entropies as well. In Ref. [21], a different approach is given and 
positive entropies are found for the phantom fluid by assuming that its temperature 
is negative defined. Babichev et al. demonstrate in Ref. [22] that the mass of a 
black hole decreases due to accretion of phantom energy. Making use of this result, in 
Ref. [23] the second law of thermodynamics {dSrotai/ dt > 0) is studied in a model of 
universe dominated by a phantom fluid in the presence of black holes (the density of 
black holes being much smaller than the density of phantom dark energy), concluding 
that it is not preserved inside the event horizon. 

Another approach applies local thermodynamics [24]. In this paper we use it 
to study a Bose-Einstein gas of self-interacting particles that produces acceleration, 
applying statistical mechanics. Thus, our model has a detailed microscopical 
description. Bose-Einstein condensates have been widely studied in cosmology. The 
better candidates for cold dark matter (axions, WIMPs, scalar fields) can be found in 
the Universe in the form of Bose-Einstein condensates [25]. Also the quantum field 
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equations describing the Bose-Einstein condensate allow to test experimentally some 
of the predictions of semiclassical quantum gravity models of general relativity [26 . 
Here we demonstrate how accelerated expansion can be obtained by considering that 
the universe contains a Bose-Einstein gas of self-interacting particles. Specifically, the 
interaction between the particles of the gas gives accelerated solutions (with p+3p < 0, 
and even phantom- like solutions with p + p < 0). It is well known in the literature 
that the presence of a self-interaction between elemental particles can originate an 
accelerated expansion of the universe [27] ■ In addition, Bose-Einstein gas is a feasible 
option to describe both dark matter and dark energy, in contrast to radiation or 
ordinary non-relativistic matter. We also demonstrate how this model can be applied 
to describe the early stage of accelerated expansion, solving the horizon problem and 
respecting the nucleosynthesis scenario. In a future paper, we will demonstrate how 
present day accelerated expansion can be achieved in our model. 

In Sec. II, we use expansions within statistical mechanics to obtain the equations 
of state of this gas. In Sec. Ill, we study the dynamics of a Friedman-Robertson- 
Walker (FRW) universe that contains the interacting Bose gas, making use of the 
results of the previous section. We obtain the parameter conditions leading to 
accelerated non-phantom and phantom solutions. In particular, we study the dynamics 
of a FRW containing the gas, and then with radiation. Then, we show how the 
horizon problem is solved for solutions with p + p < 0. We also investigate how the 
model matches the standard cosmology, respecting the nucleosynthesis scenario, and 
producing adequate fluctuations. Finally, we study the stability of the model. In Sec. 
IV, we summarize our findings. 

We use units with H = fcs = c = 1. 



2. Short-range interactive Bose gas 



The thermodynamic potential f2 of a Bose gas, described by single-particle quantum 
states with energy levels e k associated to momentum k [28], is given by 



= T^log(l-e^)/ T ), (1) 

k 

where T is the temperature, and p the chemical potential. It leads to the average 
occupation number 

fi * = (2) 
where e k contains the kinetic energy of a particle with momentum k and mass m; 
considering phase space, the number of particles associated to a given state is 

dN k = 97K-v^n k d 3 k, (3) 

where for spin-zero particles the degeneracy factor is g — 1 and V is the volume. The 
total number of particles is 

N = J dN k . (4) 

Consistency of Eq. (0} requires macroscopic occupation of the lowest state at low 
enough temperature T. In the free case, this starts at the critical temperature T c such 
that 
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or equivalently 

The chemical potential remains constant and zero for T < T c ; the latter condition 
is required for the condensate phase to be present. In addition, a non-relativistic 
description is assured if m 3> T c , or equivalently, 

A short-range two-particle interaction V(x) modifies free-particle behavior. We 
consider temperatures below the critical temperature, T < T c . Near T c , and T < T c , 
the single-particle energy is given (to first-order) by Ref. [31] 

k 2 

£k = -z h2z)on e>0 , (8) 

2m 

where the second term is the potential energy of N e> o non-condensate particles 
associated with V(x), n e> o — N e>0 /V, and 



d 3 xV(x). (9) 
For Eq. ©, the chemical potential is 



H = 2n c v' Q , (10) 

instead of the fj, = value for a free gas with T < T c . The value of T c in Eq. (|5]) and 
the constancy of \x remain valid also in this case. 

A similar constant potential term is approximated near total condensation T ~ 
[28] with n c — N c /V, N c the condensate particle number, for which 

k 2 

and 

fx = n c t>Q. (12) 

The presence of a factor 2 within Eq. (|10|) expresses the exchange effect, which 
accounts for interaction of different states, while these corrections are absent to this 
order for the condensate. We use the limits in Eq. ((SJ) for T near T c and Eq. (TlT]) 
for T near to construct an interpolated potential term. This takes into account the 
potential particle exchange between the condensate and non-condensate components, 
so that an appropriate average is v = v' [(N 2 + 2(N - N c ) 2 + 2N C (N - N c )]/N 2 , 
(N = N C + N e>0 ). 

One way to derive the energy corresponding to is using the thermodynamic 
identity 
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on on 

E = -^- T 8f + a (13) 

Another way to obtain it in the interactive case is to substitute the single-particle 
distribution into the Bose distribution in Eq. ([3]). Then, 

_ gVm 3 / 2 /" , 3/2 1 1 A 

E = J ^ e(^-M)/T-i + 2V ^ U0 (W) 

„,5^m _ T 3 /2 / + ^V 2 , (15) 



V2tt 2 Jo e 2 -l 2^ 

where the second term in E sums over pairs of interactions averaged over volume from 
Eq. ©, and Eq. (fl4l takes the thermodynamic limit. 

The pressure can be obtained from the definition p = — (dE /dV) N s and is 



P= oT7 + Z u 0T72- ( 16 ) 



2E IN 2 
SV + 6 V °V 2 

We obtain the entropy of the gas S in terms of the energy E, number of particles 
TV and volume V from the well-known thermodynamic relation 

TS = E+pV- 2V>, (17) 

where the temperature is T — (dE/dS)v,N- In our case, Eq. (flT)) reads 

'dE\ „ 5 ( „ 1 TV 2 * 

V,N 

which can be integrated to 



2 / 1 7V 2 \ S 
S = CV*\E--v Q —\ , (19) 

where C = | (128<7) 2 ^ 5 m 3 / 5 is set from the interactive-gas entropy dependence on the 
kinetic energy, which is similar to that of the non-interactive gas entropy in Ref. [28 . 

3. Cosmology of the interacting Bose gas 

We consider a flat Friedmann-Robertson- Walker (FRW) universe with line element 

ds 2 = -dt 2 + a{t) 2 [dr 2 + r 2 dn 2 ] , 

where a(t) is the scale factor, t, r and fig are the time, the radius and solid-angle 
comoving coordinates of the metric, respectively. The Einstein equations |17j read 

fl 9 =f-V = A-P. (20) 
\a J dm p 

p + 3H(p + p)=0, (21) 

where mp is the Planck mass, p and p are the total energy density and total pressure 
of the fluid conforming the universe, respectively. For a universe containing n fluids, 
the total energy density and pressure are p — Y^j=o Pj ano - P ~ S T=o Pj > wnere Pj an d 
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Pj are the energy density and pressure of each component j . If there is no interaction 
between the different components, we can generalize Eq. (|21[) to a set of n independent 
equations 

p j + 3H(p j +p j )=0. (22) 
Combining Eqs. ([20)) and (|2~Tj) , we find the acceleration parameter 

d 47T , , „ 

(p + 3p). (23) 



a 3njp 

One concludes from this equation that the universe is undergoing an accelerated 
expansion if the fluid violates the strong energy condition p + 3p > 0. Additionally, 
we can consider the dominant energy condition p + p > that assures that the speed 
of sound in the fluid be lower than c. 



P=oP+7 v 077J- ( 24 ) 



3.1. Friedmann- Roberts on- Walker metric containing the interacting Bose gas 
The pressure of the gas, from Eq. ((T9|) . is 

1 N~ 
3 r 6 V 

In an expanding FRW universe the volume changes with time as V = Vi (a(i)/aj) 3 , 
where Vi and Oj denote the volume and scale factor at an arbitrary reference instant ij , 
respectively. Assuming a FRW universe model that expands adiabatically, the above 
equation reads 

*H' + r° n K5)) 6 ' (25) 

where rii is the number density of particles (i.e., ni = N/Vi, where iV is defined in Eq. 
(Q}) at ti. Using the equation of state (|2~5|) in Eq. (|2"Tj) . we obtain the evolution of the 
energy density 

^^(5y) 5+ ^(5)) 6 ' (26) 

where p C i denotes the energy density of the condensate at the reference instant ij. 

The Bose-gas equality T = (5/3)E c /S (Ref. [35]) for the condensate energy 
contribution E c = p c V confirms that once the condition T < T c is satisfied at some 
point, it is always satisfied, as both T (Eq. ([26])) and T c (Eq. ([6])) decrease as a -2 . This 
also means the non-relativistic condition in Eq. ([7]) can be consistently maintained 
throughout. 

From Eq.(|23j). the universe acceleration containing the gas is 



a 4-7T 



a 3rap 



(27) 



It is possible to obtain accelerated solutions to the Einstein equations if the interaction 
term is negative (vo < 0), i.e., if the interaction is attractive. From this point on, we 
assume < [jjj. When the acceleration given by Eq. ([2"T]) is positive at a — (i.e., 

| In this case, the second term in the right hand side of Eq. i25l can be interpreted as an effective 
viscous pressure. The dynamics of a universe containing a viscous fluid has been widely studied and 
its presence allows the expansion to accelerate under certain conditions |32l . |33l . 
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—1vtfn\ > 3p c i), the model goes through a transitory stage of accelerated expansion 
that ends when the interaction term (which scales as a" 6 ) becomes negligible as 
compared to the condensate term (which scales as a~ 5 ). We evaluate the scale factor 
for which the accelerated stage ends (i.e., a(a/ai) = 0). This is at 

± = (28) 

Oi 3 p c i 

General Relativity (GR) requires that the total energy density be positive, i.e., 
—vonf < 2p C i. Using this condition in Eq. ([28]) , we conclude that the accelerated 
expansion stage ends at a scale factor a/ai always lower than 4/3. The transitory 
acceleration stage predicted by this model is hence more suitable to describe the 
acceleration era of the early Universe than the present accelerated era. To describe 
the present era of expansion with this model, the interaction term should evolve with 
the scale factor in a different way to become the dominant term of the expansion later 
in time. This will occur if the interaction parameter vq evolves with time, or if the 
number density of the gas particles n is not conserved with the expansion. 

We now explore the conditions under which the transitory accelerated stage 
predicted by this model can describe the acceleration of the early Universe, and in 
particular, solve the horizon problem. 



3.2. Friedmann-Robertson-Walker metric containing the interacting Bose gas and 
radiation 

To address this application of the accelerated stage, we add to the model a radiation 
term with density p r = p r i (ai/a(t)) A and pressure p r = p r /3- This term is introduced 
to fit the standard cosmology, as indeed it eventually dominates the expansion, 
allowing for nucleosynthesis to happen. 

The total energy density and pressure of the universe are 

\ 4 / \ 5 i / 



a(t)J +Z 2 V ° n *{W)) < (29) 



1 / CLj, \ 2 I 0>i \ 1 2 < ai 



P =3 Pri {W)) + 3^U)J + 2 V0nt {W)) ■ m 

The free parameters in this model are the density of the radiation p r i, the energy 
density of the condensate p C i, and the interaction term Von?. The parameters are 
constrained by the condition —vonf < 2(p r i + p C i), as the total energy density entering 
the Einstein equations is required to be positive, i.e., p(ai) > 0. The reference term 
ai can be interpreted as the scale factor at which our model starts to be valid. 
The acceleration reads 



a 4-7T 



a 3m|, 



• 1 TTTTT.) +ipc \w)) +2v ^\W) 



(31) 



Similarly to the model in subsection 3.1, it is possible to obtain accelerated solutions 
to the Einstein equations if Vo < 0. The acceleration eventually becomes negative, as 
the interaction term scales as (ai/a) 6 , while the contributions of the radiation and the 
condensate scale as (ai/a) 4 and (ai/a) 5 , respectively. When the latter terms dominate, 
the expansion becomes decelerated. 
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Figure 1. Left panel: Energy-density evolution from the reference scale factor 
di on, for different choices of the free parameters. Right panel: Acceleration 
evolution with the scale factor from the reference scale ai on, for different choices of 
the free parameters. Every line is normalized to the value of the density at a = ai, 
Pi = Pri + Pci + (^0ii)/2. The three different regimes are shown: i. Decelerated 
universe with parameters p r i = 0.5pi, p c i = 1.5pi and (von?)/2 = —pi (solid 
line); ii. Accelerated solution with parameters p r i = 0.5pi, p c i = 5pi and 
(urjn?)/2 = — 4.5pi (dotted line); iii. Super-exponential accelerated solution with 
parameters p r i = 3pi, p c i = 5pi and (von^)/2 = —7pi (dashed line). We assume 
for this and the next plots, 87r/(3rrip) = 1. 



The left and right panels of Fig. [T] show the evolution of the total density and 
acceleration, respectively, for three possible types of solutions, which are characterized 
by the acceleration at a,. These are: 

i. Decelerated solutions 

The acceleration given by Eq. (J3J) is negative (or zero) at the scale factor a = a L 
when p(ai)+3p(ai) > 0, i.e., —v^nf < (3/2)p C i+p r i. Then, the universe maintains 
a decelerated expansion from a — ai on. 

ii. Non-phantom accelerated solutions 

These solutions are obtained when the free parameters fulfil the conditions 
(3/2)p c ,; + p r i < —VQnf < (5/3)p C i + (4/3)p r ,, where the lower limit is given 
by the positiveness of the acceleration in Eq. (|31[) at the scale factor a = a, , and 
the upper limit by the condition p + p > 0. 

iii. Phantom-accelerated solutions 

The acceleration given by Eq. (|3"Tj) at a; is positive, and the condition p + p > 
is violated at the scale a — ai, if — vqu\ > (5/3)p C j + (4/3)p r i. In this case, the 
total density grows with the scale factor as dp/ da = —3(p + p)/a > 0. First, the 
universe goes through a super-exponential acceleration phase similar to that in 
the phantom-energy models. When the interaction term is diluted for large values 
of a, the condition p + p < is no longer valid and the energy density becomes a 
decreasing function of a. Consequently, the energy density presents a maximum 
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. The scale factor a* for which p(a„) = p. 
a* — Qvonf 



max 



IS 



(32) 



a i 5p cl + v / 25p 2 cl - A8v nfp, 



The plausibility and stability of phantom models has been widely studied [34] . 
The big-rip problem (infinite expansion of the universe in a finite-time span) in 
some phantom-energy models that explain the current Universe acceleration [2] 
is not present in our model because the phantom-accelerated stage is transitory 
and ends before the scale factor diverges. 

The scale factor a at which the accelerated expansion stage ends in this scenario is 
always lower than \[2ai (see the Appendix A for details). 

3.3. Solving the horizon problem through the super- exponential expansion regime 

Next, after briefly reviewing the horizon problem, we obtain the parameter conditions 
which solve horizon problem, and are compatible with the nucleosynthesis scenario. 
Also, the spectrum and evolution of the initial density perturbations is addressed later. 

Cosmic microwave background (CMB) photons from different regions of the sky 
are in thermal equilibrium at temperatures with relative fluctuations up to 1CP 5 . The 
most natural explanation for this phenomenon is that photons from different regions 
had been thermalized by being in causal contact previously. However, the dynamics 
of a FRW universe containing only radiation and/or non-relativistic matter (dust) 
cannot describe such a situation as there was no time for those regions to interact 
before the photons were emitted. In other words, the distance light could travel from 
the instant the photons were emitted t e until the photons were released at the instant 
tdec is much smaller than the horizon distance today (to), i.e., 



where ao represents scale factor today. This problem is known in the literature as the 
horizon problem [17], [18] , 

An alternative description of this problem considers the Hubble radius 1/H, where 
H is defined by Eq. (|20|) . The Hubble radius is closely related to the horizon [35;. 
A point in space is no longer in causal contact with an observer when its distance 
to it, A, becomes larger than 1/H. In radiation and dust dominated universes, 1/H 
increases with a slope with respect to a larger than one. This implies that a point A 
in space for which A^ = 1/Hq (in causal contact with the observer for the first time 
today to) was never before in causal contact with the observer (i.e., A^i > 1/H(a) for 
any a < ao). 

In both the accelerated and phantom-accelerated cases of the FRW model 
containing radiation and the Bose gas exposed in the previous subsection, there is 
an initial period for which different length scales A are in casual contact, leaving 
the Hubble radius and reentering it later in time, when the universe is no longer 
accelerated. This process is illustrated in Fig. [2] For the accelerated region, 
d(l / H) / (da) is lower than 1 but positive until the expansion becomes decelerated. For 
the super-exponential accelerated region, d(l/H)/(da) is negative until the expansion 
becomes decelerated. This allows for a large range of length scales to leave the causal 
contact region and to reenter it later. Finally, we also note that for the decelerated 
region, this phenomenon is not manifested, as d(l/H)/(da) is always larger than one 




(33) 
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Figure 2. Two different length scales that evolve with scale factor as A = 
\j(a(t)/ai) (black solid lines). The Hubble radius, l/H, is plotted for different 
choices of the parameters: 1. Decelerated universe with parameters p r i = 0.5pi, 
p c i = 1.5pi and (von?)/2 = —pi (solid line); 2. Accelerated solution with 
parameters p r i = 0.5pi, p c i = 5pi and (von^)/2 = — 4.5pi (dotted line); 3. 
Super-exponential accelerated solution with parameters p r i = 3pi, pd = 5pi 
and (vonf)/2 = —7pi (dashed line). For this plot, pi = p r i + p c i + (uon?)/2. 



(the particle horizon tends to zero when a decreases, as for the radiation or the non- 
relativistic matter FRW models). 

It is possible to solve the horizon problem if we consider phantom-accelerated 
values for the free parameters of our model; specifically, we choose the scale factor Oj 
to be small enough so that the Hubble radius is big enough. In this scenario, choosing 
1/Hi 3> l/H (a*) (a*/aj defined by Eq. (|32)) ). the distance that light can travel before 
the photons were released can be large enough to solve the horizon problem because 

da' f to dt' 

» / ~~ 777T • ( 34 ) 



a' 2 H(a') J tdcc a(t') 

Then, we ensure that there be wide range of scales that leave the region of causal 
contact. Fig. [3] illustrates an example of a phantom regime for which 1/Hi = 
10 30 /if(a*). This implies that the energy density increases from pi until it reaches its 
maximum p max sa 10 6 V;. 

We evaluate the range of the free parameters that solve the horizon problem. It 
is plausible to assume that p max cannot be larger than the Planck scale, mf> (GR is 
valid for p < m,p). We first define the new set of parameters 




x * _ Pn I \ x* = Pci ( — ) 1 •") 

a * / Pmax \ a * 

i 2 / \ 6 

1 u n| / a, x 



= 1? J- (36) 

& l J max \ •** / 

X* has been defined to be positive when the interaction is attractive. These 

parameters are the contribution of each energy component at a* to p ma x- We define 
also 

X(a) = (37) 

Pmax 

as the total energy density reescaled to p max - Its evolution in terms of the scale factor 
can be expressed as 

*<.) (38) 

The rescaled parameters are related because X(a*) = 1, that is, 

X*+X*-X* = l. (39) 



Accelerated expansion of a universe containing a self-interacting BE gas 



12 




Figure 4. The energy density as a function of the scale factor. The shaded zone 
represents the shapes of the energy density depending on the value free parameter 
X*. The contours of the shaded zone are the limit values of the free parameter: 
when the condensate term is null (X* = 2); and when the radiation term is null 
{X* = 5). 



In addition, the derivative of the total density respect to a evaluated at a* is equal to 
zero. This condition reads 

- AX; - hX* + 6X* = 0. (40) 

Combining conditions (f3!))) and (|4"P|) . we relate the parameters in the form 

X* =5-X*, X* = 2X* - 4, (41) 

and, then, the only free parameter to consider is the interaction X*. Given Eq. (l4Tj) 



and because X* > 0, X* < 5. Additionally, X* > 0, which implies that X* > 2. The 
range of parameters that can solve the horizon problem in our model are then given 
by Eq. (j4Tj) with 2 < X* < 5. Fig. |4] shows the energy density as a function of the 
scale factor, for this range of parameters. The energy density reaches the null value 
at the scale 

a** = X* 

a* X*-2+ ^XjTZ 

and the horizon problem is solved, and Eq. [33] is satisfied, because the energy density 
increases from p(a) w at a scale factor a > a**. Note that if we consider a** as the 
initial scale, l/JJ(a**) diverges, and so does the horizon. This divergence is also found 
in the inflationary de-Sitter models, for which the horizon diverges when a(t e ) — > 0. 
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3.4- Fluctuations 

One of the main predictions of inflation is the presence of scale-invariant fluctuations 
which form the seed of later large-scale structure formation. In the conventional 
description for a nearly de Sitter universe, the constancy of the fluctuations amplitude 
is determined by the constant horizon size (which in turn is set by the Hawking 
temperature Th = H/(2tv)); this occurs at horizon scale k = Ha, where k is the 
Fourier-transformed comoving coordinate. 

We have shown that the universe can be initially causally connected in our model; 
we shall also argue that the properties of a Bose gas allow for a scale-invariant 
perturbations regime: we show that thermal fluctuations originating in the non- 
condensate contribution have the appropriate properties for the standard cosmology. 

During the phantom-accelerated regime, the expansion is dominated first by the 
interaction term in Eq. (|29[) and the non-condensate contribution; the latter then 
dominates as a grows (see Fig. 1). In general, the mean square density fluctuation, 
for given volume (and therefore, constant particle number or entropy), is given by |29j 



((E Ef) = (*^) = T 2 iVc„ (43) 

where c v is the specific heat per particle (constant in our case). 

While this calculation represents a useful benchmark for later comparison, 
we actually need the contribution of the fluctuations at given mode k. This 
may be obtained from the approximate large-distance non-condensate wave-function 
component 5$k [30] 



m= 2^ PmV £ a&v.l***! 2 ' ( 44 ) 

where p m is the non-condensate mass density, and we used the physical kph ya \ this 
term contributes to the kinetic term in Eq. (f2U)) . Then, one obtains the expresion for 
the statistical average 

W) = vT^— (45) 

V Pm^ phys 

By considering the phase-space factor Vd 3 k p h ys — d 3 k, one finds the k contribution 
to the total fluctuation in Eq. (|43|) is proportional to T. While T is a-dependent, 
the ratio in the rms density contrast is not, nor is Crms = Sp rms /(p + p), as one can 
use the relation p c — (2/3)p c for the non-condensate part. We hence obtain that the 
gauge-invariant quantity contribution £ rms is constant beyond the horizon. This leads 
to gaussian scale-invariant fluctuations, since their origin is thermodynamic, and they 
all have the same fc-power dependence for their amplitude. 

Once the fluctuation mode k crosses outside the Hubble radius, it evolves as 
a classical perturbation |17j : the fluctuation evolution of the universe in our model 
becomes equivalent to the standard universe. Thus, we find that it is possible to obtain 
scale-invariant fluctuations, assuming the non-condensate fluctuations dominate, even 
as the horizon size 1/H is not constant, but decreases with a (or equivalently with t). 
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3.5. Nucleosynthesis 

We find the condition under which our model solves the horizon problem consistently 
with standard cosmological requirements, and specifically with nucleosynthesis |17j . 
Nucleosynthesis requires that the Universe be dominated by radiation at times of 
around Is in order for the light elements to be produced. 

We consider the scale factor a eq at which the radiation and Bose-gas densities are 
equal, i. e., 

\ 4 

a* 



eq 



x; — =x;\ — \ -x; — . (46) 



From the scale a eq on, the radiation term will dominate the expansion of the universe. 
We evaluate the time span t eq which corresponds to this scale factor. If t eq < Is then 
the nucleosynthesis scenario will be preserved by the model. The time span t eq is 
found by integration of Eq. ([20"]) . 




-a d ta_ 



(47) 



If there is no radiation in our model, the parameters read X* = 0, X* = 6, 
X* = 5. In this particular case, t eq tends to infinity. These parameters, although they 
solve the horizon problem, are in contradiction with the nucleosynthesis scenario. 

Depending on the value of p ma x it is possible to numerically evaluate how much 
radiation density is initially needed (i.e., how large X* must be) to allow for t eq < Is. 
This minimal radiation amount decreases its value as p max increases, as can be seen 
from Eq. (|i?)t . 

Assuming p max = m% and that X* w 1(T 17 , I*w6-2x 1CT 17 , X* « 5- 1CT 17 , 
then Eq. (|4T|) reads t eq rj lO 43 ??!^ 1 rj 0.1s. In this limit, nucleosynthesis will be 
allowed for any X* > 10~ 17 . If we consider that p max is lower than rrip, the minimal 
amount of radiation needed (evaluated from Eq. (I4T[) with t eq = Is) is larger than 

io- 17 . 

Consideration of the Bose-gas particles' mass requires the contribution mn in 
the energy density. This term scales as a -3 , and eventually becomes the dominant 
contribution. To ensure that this term does not dominate the expansion before 
nucleosynthesis, we have to impose the condition mn, <C pri- The latter condition, 
together with the non-relativistic gas condition given by equation ([7]). implies p„ -C 
pri, (vonf)/2. In terms of the rescaled parameters, this is X* s» 0, X* sa 2 and X* ps 3. 

We can conclude that nucleosynthesis does not bound the free parameters of the 
model strongly. If the radiation term does not differ from the other two terms for more 
than 10 orders of magnitude, nucleosynthesis will occur. In such a scenario, our model 
starts with an accelerated era that explains the horizon problem, later evolves in a 
radiation dominated universe in which light elements can be created and, eventually, 
the expansion is dominated by non-relativistic matter. In other words, our model of 
universe behaves as the standard cosmology, besides the initial stage of expansion. 
As we mentioned before, late time acceleration can be also obtained in our model if 
we consider that the number of particles of the gas or the interaction term are not 
constant. This problem will be addressed in a future paper. 
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3.6. Stability of the interacting Bose gas in the cosmological frame 

A confined Bose-gas condensate with an attractive interaction is stable under certain 
conditions [35], unlike the case of a Bose condensate with an attractive interaction 
within an infinite volume, which always collapses. When the gas' particles arc confined 
in a given volume (e.g., by a harmonic trap), finite-size effects due to the confining 
potential generate an additional kinetic-energy contribution that compensates the 
interaction energy, avoiding the collapse, as long as the confined particle number 
in the potential is lower than a critical value N cr . This number can be computed by 
equating the exterior potential plus kinetic energy terms with the interaction term 
[36] . which translates to 

N cr < 4tt^, (48) 
m\v \ 

where ano is the physical trap characteristic length scale, and m is the particle mass. 

We approach the collapse problem in our model applying dimensional analysis. 
In the cosmological frame, we assign a% to the Hubble radius 1/H, as the latter 
represents the region of space in causal contact. In our case, aho grows with time, 
as the volume expands [37]. Thus, if the Hubble volume (1/H) 3 contains less than 
N cr particles, the collapse is avoided. So is the case for any smaller volume with a 
characteristic length ahoi lower than l/H, as such a volume will contain less than N cr 
particles. 

It is sufficient to apply the above condition for initial values of the parameters in 
our model, as the interaction energy term (proportional to n 2 ) dilutes faster with the 
expansion than the condensate term in Eq. ff29|) . and eventually becomes negligible. 
With the association of N cr to the ground-state (or condensate) particle density n c , 
Eq. dH]) reads 

2 + 2(47r) 2 nJ 

where we used a^ G ~ l/-ff; = {[8tt / (3mp)](p ri + p C i + uonf/2)} _1 / 2 , and m is the total 
initial particle density. Using the non-relativistic condition for the mass in Eq. (|7|). 

1/3 

m ^> rij , and having in mind that m < mp , it is convenient to explore two extreme 
cases 

i. Hi ^> n c 

Only a small amount of the total particles are in the ground state. As the 

temperature of the gas is T = T c (^^) 2/3 [28], we have T ~ T c . If m is 
such that the second term in the denominator in Eq. (|49[) satisfies 

then Eq. ([4*$)) reads n^|v | ^ 2(p ri + p ci ), which is consistent with the condition 
required for the total energy density entering the Einstein equation to be 
positive. In this case, it is possible to obtain phantom-accelerated, non-phantom 
accelerated and decelerated solutions that do not collapse. 
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All the particles of the gas are in the ground state (i.e., the temperature of the 
gas is T = 0,) and consequently, 

3m p mn c 3 mpm 

2(47r) 2 nf = 2(4^)2 > 1 

2 + 2(47r) 2 n^ 

In this case nf \vq\ -C + p C i. The condition of p > is compatible with 
the latter condition. But conditions for the accelerated and phantom kind 
solutions in Section 3.2 ((3/2) ( o ci + p ri < —v nf < (5/3)p c , + (4/3)/?h and 
—Vorq > (5/3)p C i + (4/3)p r j, respectively) are in contradiction with the collapse 
condition in this case. Consequently, models for which m — n c and that do not 
present collapse are suitable to be used in cosmology (as p > 0) but lead always 
to a decelerated expansion of the universe. 

The above conditions are further relaxed in a cosmological context, as an attractive 
interaction ultimately accelerates the expansion, and temperature effects will 
contribute to counter the attraction. To summarize this subsection, stable solutions 
of the Bose-Einstein gas are possible, while accelerated ones favor the case m 3> n c . 

4. Conclusions 

A Bose-Einstein gas with self-interacting particles leads to accelerated expansion of 
the universe. The interaction between the particles of the gas gives naturally solutions 
of the type p + 3p < 0, if the interaction is attractive. Phantom-accelerated solutions 
with p + p < can be also obtained, if we consider a model of universe that contains 
the Bose gas plus radiation and in which the interaction term initially dominates the 
expansion. 

As such a term, responsible for the accelerated behavior, decreases faster with the 
expansion than the condensate and radiation terms, the accelerated stage eventually 
ends. This fact makes the model more suitable to describe the acceleration of the 
early Universe than the present accelerated era. 

Our model solves the horizon problem when we consider phantom-accelerated 
solutions. In addition, nucleosynthesis is feasible by starting with a large-enough 
radiation term. 
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Appendix A. 

In this Appendix we calculate a higher limit for the scale factor a for which the 
accelerated stage of the FRW model of interacting gas plus radiation of Sec. II ends. 
Assuming that the acceleration of the expansion given by Eq. (l3"Tj) is positive at the 
scale ai, we thus evaluate the scale factor a for which a(a) = 0. This is at 
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By definition, p ri > 0. We define the new variables 

y=^>0, z = - v jdll>0, (A.2) 

Pri Pri 

and express Eq. (IA.1[) in terms of them 



a _ -3y + y/9y 2 + 32z 



(A.3) 



GR requires the total energy density to be positive (—vgnf < p r i + p c i), which reads 
in terms of the new variables z < 1 + y. Making use of the latter condition and of Eq. 
(|A.3|) . we get 



< _3y + W + 32( y+l)_ 



4 



The right-hand side term of inequality (|A.4j) is a decreasing function of y. Its maximum 
value in the limit y — > is Consequently, we conclude that the acceleration stage 
in the radiation plus interacting Bose gas scenario ends at a/ai always lower than y/2. 
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